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SUMMARY 

Superposition techniques are used to calculate the rate of heat 
transfer from a flat plate to a turbulent incompressible boundary layer 
for several cases of variable surface temperature. The predictions of 
a number of these calculations are compared with experimental heat- 
transfer rates, and good agreement is obtained. A simple computing 
procedure for determining the heat-transfer rates from surfaces with 
arbitrary wall- temperature distributions is presented and illustrated 
by two examples. The inverse problem of determining the temperature 
distribution from an arbitrarily prescribed heat flux is also treated 
both experimentally and analytically. 


INTRODUCTION 

report is the third in a series of four covering a three— year 
investigation of heat transfer to the turbulent incompressible boundary 
layer with arbitrary wall -temperature variation (ref. l).l In the first 
report the experimental apparatus is described, and the results of ex- 
periments with constant wall temperature are given (ref. 3). The second 
report presents the results of experiments and analyses for a step wall- 
temperature distribution (ref. 4). In the present report the step- 
function analysis is used as a basis for predicting heat-transfer rates 
for several cases of variable wall temperature, and the predictions are 
compared with experimental data. A simple method for handling arbitrary- 
wall- temperature problems analytically is also presented. The fourth 
report presents an analysis of the effect of transition point on heat 
transfer in the turbulent boundary layer and compares the results with 
experiments (ref. 5). 


■kPhis nonis othermal heat-transfer work is summarized briefly in 
reference 2 . 
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Because of the linearity of the boundary- layer energy equation, the 
heat -transfer characteristics for a step wall- temperature distribution 
may be superimposed to determine heat-transfer rates for arbitrary- wall- 
temperature situations. This fact has been pointed out by Klein and 
Tribus (ref. 6) and by Rubesin (ref. 7). The superposition results in 
integrals that must be evaluated at each point where the heat-transfer 
rate is desired. The integrations are easily performed if the wall- 
temperature variation is relatively simple, such as linear, parabolic, 
and so forth, but numerical methods are required for more complex cases . 
One approximate method for treating variable-wall-temperature problems 
is to express the wall-temperature distribution by a finite number of 
steps. However, this results in infinite heat-transfer rates at the 
steps and thus does not yield meaningful results in the region near the 
discontinuities. A better method, short of exact integration, is to ap- 
proximate the temperature distribution by a finite number of linear seg- 
ments or "ramps." Since this approximation is continuous, the heat- 
transfer rates are everywhere finite, and good results can be obtained 
with relatively few ramp segments. This technique is described in the 
present report, and the results compare favorably with experiments. 

The inverse problem, wherein the heat-transfer rates are specified 
and the temperature distribution is to be determined, can also be handled 
by superposition of steps or ramps. This problem is also discussed in 
the present report. 

This investigation was carried out at Stanford University under the 
sponsorship and with the financial assistance of the National Advisory 
Committee for Aeronautics. 


SYMBOLS 


A(a/x) function for prescribed temperature problem, 
(l0/9)B r (8/9,10/9) - [l - (a/x)] 

a location of ramp, ft 

b height of step, °F or Btu/(hr)(sq ft) 

c p specific heat at constant pressure, Btu/(lb)(°F) 

D(a/x) function for prescribed heat-flux problem, 

A a\ B r (l/9,20/9) - (a/x)B r (l/9,10/9) 

V x) f(i/9)r(a/9) 

G free -stream mass velocity, pu^, lb/(hr)(sq ft) 
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gfc ,x) kernel function for prescribed heat flux 
H( Z/x) function for prescribed heat-flux problem, 1 - 

h local heat -transfer coefficient, q"/At, Btu/(hr)(sq ft)(°F) 

h(£,x) local heat -transfer coefficient at x due to a step temperature 
at £ 

, k thermal conductivity, Btu/(hr) (ft ) (°F) 

l location of step, ft 

m slope of ramp, °F/ft or Btu/(hr)(cu ft) 

Pr Prandtl number, pCp/k 

d q" heat flux, Btu/(hr)(sq ft) 

■3 

3 q* effective q", Btu/(hr)(sq ft) 

i 

5 - Re x Reynolds number based on x, Gx/p 

r 1 - (a/x) 9 / 10 or 1 - (z/x) 9 / 10 

S(z/x) function for prescribed temperature problem, 

9/101-1/9 

.1 - (l/x) J - 1 

St local Stanton number, h/Gcp 

St jj local Stanton number for constant heat input 

St^ local Stanton number for isothermal plate 

T w absolute plate temperature, °R 

absolute free- stream temperature, °R 
t temperature, °F 

At t w - t„, °F 

At* 


B r (l/9,10/9) 

r(i/9)r(8/9) 


effective At, °F 

mean temperature of heated strip, °F 
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V - t„, °F 

plate temperature, °F 

free -stream temperature, °F 

free- stream velocity, ft/sec 

load distribution 

distance from leading edge, ft 

distance from plate, ft 

variable of integration 


B(a,b) 

beta function, 

f ^-Ri- 

0 

z) b_1 dz 

B r (a,b) 

beta function, 

A 

f ^-Ri- 

D 

z) 13 " 1 dz 

r(a) 

gamma function, 

v* 00 

J e -z z a-l 

i — i 

i 

dS 

ii 
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M- 

£ 

p 

<P 
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viscosity, lb/(hr)(ft) 
variable of integration 
density, lb/cu ft 
beam deflection, in. 
beam influence coefficient, in. /lb 

QUALITATIVE EFFECTS OF VARIABLE WALL TEMPERATURE 

Until quite recently, practically all Investigations of boundarv 

ature Mth^hRh treated . the s P ecial case of constant surface temper- 
ature. Although the assumption of constant temperature greatly simnli- 

les analysis, many important systems involve heat transfer from non- 

iso ermal surfaces, and failure to consider the effects of the noniso- 

rat^T Recent effoVT t0 err ° rS “ calculat ed heat-transfer 

rates Recent efforts have provided methods for treating the noniso- 

thermal problem, and the calculation of heat-transfer rates is now es- 
sentially mathematical. 



4994 


5 



Before investigating the details of calculation of heat-transfer 
rates from nonisothermal surfaces, it is desirable to obtain a qualita- 
tive understanding of the way nonisothermality might be expected to in- 
fluence the heat transfer. This can be best achieved by examining the 
istory of the boundary layer and considering the qualitative effect 
of events upstream on the temperature profile in the boundary layer. 

ile these remarks apply to a heated plate, they may readily be extended 


If the plate is at constant temperature, the temperature profiles 
at different points on the plate will have the same general shape and 
will appear as shown in sketch (a). The slope at the wall will be pro 
portional to the heat-transfer rate. " 



If the plate temperature increases in the flow direction, the tempera- 
ture profiles will tend to be more drawn out, as indicated in sketch (b). 
The gradient at the wall will therefore be steeper, and the heat-transfer 
coefficient will be greater than if the plate were at constant 
temperature . 



(b) 
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If, on the other hand the^U -P-t^ec^aeee 

tion, the profiles w therefore he less, and. the heat- 

S S -fS-e were at constant temperature. 



In fact, 
enough , 
would, he 
ference . 


t is conceivable that if the 

e profile could be so distort ® f ^ te 0 ver-all temperature dif- 
^IrcTS; te^erature pr^e^uld appear as follows: 



(d) 


Moreover, decreasing the -11 te-perature even g^^g^^which " 

complete reversal of the proil lth positive over-all tem- 

case the heat transfer would be negative, win 

perature difference. 
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It is also possible that decreasing the temperature difference even 
faster could lead to finite negative heat transfer with a zero tempera- 
ture difference: 



It is evident from these considerations that the nonisothermality 
may have a profound influence on the shape of the temperature profiles 
in the boundary layer and consequently may greatly influence the heat 
transfer. Moreover, it is clear that the ,f history M of the boundary layer 
is of great importance. It should also be noted that the heat-transfer 
rate need not always be in phase 1 with the wall -temperature variation. 

In general , however , the following is true (if t w > t ): 

(1) A decreasing temperature difference leads to heat -transfer rates 
that are lower than those for an isothermal plate. 

(2) An increasing temperature difference leads to heat-transfer 
rates that are higher than those for an isothermal plate. 

With the effects of thermal history well in mind, the calculation of 
heat transfer from nonisothermal surfaces becomes purely a mathematical 
matter. 
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THEORY AND ANALYSIS 
Review of Theory 

The methods for determining heat transfer from nonisothermal sur- 
faces are similar to the methods used in determining the deflection of 
beams subject to arbitrary load distributions. The energy equation of 
the boundary layer is linear in the fluid temperature if fluid properties 
are assumed to be constant This allows superposition techniques to be 
employed. Rubesin (ref. 7) has shown that the heat-transfer rate for an 
arbitrary wall- temperature variation can be determined by superimposing 
a number of. step wall- temperature distributions/' so that summation of 
the steps yields the actual variable temperature distribution, and the 
heat-transfer rate at any point is equal to the sum of the heat-transfer 
rates attributable to all "steps" upstream of the point in question. 

This idea of superposition is illustrated in figure 1, where the super- 
position of temperature steps is compared with the superposition of point 
loads used in beam- deflect ion problems. It is evident that a satisfac- 
tory solution for a step wall -temperature distribution is required before 
any attempt can be made to handle the variable-temperature problem, and 
such a solution for turbulent incompressible flow over a flat plate with 
a step wall-temperature distribution is presented in reference 4. 

It is convenient at this point to introduce some new notation. 

Since the temperature distribution along the heated plate may be thought 
of as a function of the distance from the leading edge x, one may write 

At = At(x) 

On the other hand, the temperature distribution can always be represented 
by some algebraic expression involving x, and this expression may be a 
function of several parameters. Then one might prefer to denote the tem- 
perature difference more completely by 

At - At (a-^a^a^, . . *J a n j x ) 

where the a ± are the parameters in the functional description of the 
temperature difference. They might be the coefficients in a power series 
expansion, the locations of discontinuities, or some other parameters. 
Hereafter the parameters will always be listed before the important vari- 
able, which in this case is x. If only one symbol appears in the paren- 
theses, it will refer to the important variable; and the fact that the 
function depends on the parameters aj_ to a^ is to be understood. 

Thus, for a step wall-temperature distribution, the heat-transfer coef- 
ficient may be written as h(Z;x), where l is the location of the dis- 
continuity, or merely as h(x), where it is to be understood from the 
context that the coefficient refers to the step distribution case. 
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It is also convenient to compare the local rate of heat transfer 
from a plate having variable wall temperature with the rate that will 
occur if the surface temperature is constant at its local value from the 
start of the plate. For this purpose it is convenient to denote the 
local Stanton number for an isothermal plate by St^. The isothermal 
Stanton number is a function of x only and does not depend on any of 
the parameters characterizing the nonisothermal problem. The notation 
Strp(x) will sometimes be used to emphasize this point. The isothermal 
Stanton number may be taken from any suitable expression for the Stanton 
number for heat transfer from a plate at constant temperature. For ex- 
ample, for heat transfer to a gas in the turbulent incompressible bound- 
ary layer, the isothermal Stanton number may be determined from 

St T (x) = 0.0296Pr'°- 4 Re"°- 2 f^ (l) 

Reference 3 shows that this relation is satisfactory for gases in the 
Reynolds number range 10 5 < Re x < 10^. In using equation (l), the 

fluid properties in the Stanton, Prandtl, and Reynolds numbers are eval- 
uated at the free- stream static temperature; the factor (T w /T 00 )“ (9, 4 
provides the temperature-dependent fluid-property correction. 

The step wall-temperature distribution case, which is the entire 
basis for superposition in arbitrary- wall-temperature problems, is 
treated in reference 1. The step wall- temperature distribution may be 
written as 


At = At(Z;x) =0 x < Z 

At = At(Z;x) = AtQ x > Z 


Reference 4 shows that the corresponding heat-transfer rate may be repre- 
sented by 


St ( Z ; x 
Strp(x) 



x > Z 


( 2 ) 


The heat -transfer coefficient for the step case is therefore 

9/ldl" 1 / 9 


r 


h(Z;x) = Gc-pStqi(x) 


-(D 


(2a) 


Following the methods of Klein and Tribus (ref. 6), one can superimpose 
an infinite number of small steps. This results in heat-transfer rates 


from the nonisothermal surface given by the following integral 
expression: 


1 


ii 



h(£jx) dt w (£) 


( 3 ) 


Here the kernel function h(C,x) is, from equation (2a), 


h(£:x) = Gc Stm(x) 
P 


r 

l 

_ 



( 4 ) 


Note that the terms GCpStip (x ) may all be brought outside the integral 
in equation (3), since the integration is performed over the variable £ . 
It should be noted that the integral of (3) must be taken in the 
"Stieltjes" sense (ref. 8) rather than in the ordinary "Riemann" or 
"area" sense. This must be done because the prescribed wall temperature 
may have a finite discontinuity, so that dt w is undefined at some 
point. The Stieltjes integral may, however, be expressed as the sum of 
an ordinary or Reimann integral and a term that accounts for the effect 
of the finite discontinuities. The integral of equation (3) may be 
written as (ref. 8) 



h(£jx) dt w (£) 
(Stieltjes) 



dt-.CC ) 

h(Cjx) — pp — d£ 


d£ 


(Ri 


lemann 


N 


+ d( 2 n ;x) t w (Z n ) 
n=l 



( 5 ) 


where [t (l*) - t w (Z“)J denotes the temperature rise across the n ^ 
discontinuity in the wall temperature. The use of equation (3) will be 
illustrated later by several examples. 


Equation (3) is useful if the wall-temperature distribution is pre- 
scribed and the heat flux is to be calculated. An equally important 
problem is the case in which the heat flux is given and the wall temper- 
ature is to be found. Again following Klein and Tribus (ref. 6), the 
wall temperature may be determined from 


ej? uticjis. 4yy4: 
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where 


gfejx) 


(9/10) (l/x) 

r(l/9)r(8/9)GCpSt T (x) 



(7) 


Because of the nature of the integrand of equation (6), the integration 
can always be performed in the usual "Riemann” sense. The use of equa- 
tion (6) in problems where the heat flux is prescribed will be illus- 
trated later by examples. 


Functions of Interest 

In the solution of various arbitrary-wall- temperature problems , 
integrals are frequently encountered that cannot be integrated in closed 
f orm, but they are well known and may be determined from tables . Of 
particular interest here is the beta function, defined by (ref. 9) 

B(a,b) = / z a - 1 (l - z) b_1 dz (8a) 

so 

The beta function is a function of two arguments, a and b, and is not 
tabulated as such. It is, however, related to gamma functions by the 
relation (ref. 9) 

B(a,b) = (8b) 

The gamma function is another integral but is a function of one argument 
only, and tables of this function are available (e.g., ref. 10). Note 
that B(a,b) = B(b,a). 

More properly, the beta function defined previously is referred to 
as the complete beta function. The incomplete beta function is also of 
interest, and it is defined by 

B r (a,b) = f z a " 1 (l - dz (8c) 

Jo 

Note that the complete beta function is B 1 (a,b). The incomplete beta 
function is tabulated, but not for arguments of interest in nonisothermal 
boundary- layer calculations. A number of incomplete beta functions of 
interest in turbulent heat -transfer analysis have been determined as a 
part of the current work. These functions are given in table I and 
are plotted in figure 2. In addition, the reader is referred to refer- 
ence 11 for similar tabulations of incomplete beta functions for both 
laminar and turbulent flow over a flat plate. In calculating the 
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integrals use was made of the symmetry property of the incomplete beta 
function (ref. 12) 

- B-^(a,b) - B-^_^(b,a) (8d) 

and the relation of the incomplete beta function to the hypergeometric 
function (ref. 12) 


a 

B r (a,b) = — F(a,l-b;l+a;r) (8e) 

where the hyper geometric function is a well-known series. An IBM 650 
computer was used to calculate the hypergeometric functions and the in- 
complete beta functions. 


Some Special Nonisothermal Heat-Transfer Calculations 

To illustrate the methods of superposition previously discussed, 
several examples of variable wall temperature and variable heat flux" 
have been worked out in detail. These examples illustrate not only the 
methods of calculation, but also several interesting aspects of noniso- 
thermal heat transfer. The results of these calculations are summarized 
in table II. 


Constant heat input . - The temperature distribution along a flat 
plate at constant heat input may be determined from the theory of 
variable-wall-temperature heat transfer. The temperature distribution 
is, from equations (6) and (7), 


(9/lO)q5 

W “ rU/9jr(8/9)Gc p St T (x) 
Setting z = 1 - (l/x) 9 / 10 reduces 

u 

Ax , ^ !0 

AtW ~ r(,l/9jr(8/9)Gc p St T 




(9a) 


(9b) 


The integral of (9b) is recognized as a complete beta function, and 
therefore 


, , <$>(1/9,10/9) q 0 

r(l/9)r(8/9)Gc p St T U) ~ °- 9587 GcpStrp dOa) 
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Putting (lOa) into dimensionless form gives 


St 

St 


= 1.043 


(10b) 


H ^ e notes the local Stanton number for constant heat input 


where St 

-•* x^x uuiibuano neap input 

heat S™t a i eS i ^ the t Cal heat - transfer coefficient for constant 
heat input is 4.3 percent greater than that for constant wall tempera- 
ture at any given location on the plate. P 


gamp wall temperature . - Consider the 
tribution as shown in sketch (g): 


'ramp' 1 wall-temperature dis- 



The wall temperature is given by At(x) = mx, where m is the 
the wall temperature, dAt/dx. Substituting in equation ( 3 ), 


slope of 


q"(x) = Gc p St T (x)m 1 - g) 9 / 10 ' 


, 1 - 1/9 


d£ 


(11) 


By letting z - 1 - (5/x) 9 /l°, the integral of (ll) 
beta function and leads to 


reduces to a complete 


q"(x) - 0CpStQi(x)mx (l0/9 )B( 8/9, 10/9) 


= 1 . 134mxGCpStip (x ) 


( 12 ) 


Since mx is simply the local temperature difference 
be written as 9 


equation (l2) may 


St 

st T 


1.134 


(12a) 
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Hence, the heat-transfer rate is 13.4 percent greater than would be pre- 
dicted from the equation for a const ant -temperature plate. 

Step-ramp wall temperature . - Consider the wall- temperature distri- 
bution indicated in sketch (h): 



The heat-transfer rate for this wall- temperature distribution may be de- 
termined by superimposing the heat-transfer rates for a step at the lead- 
ing edge and a ramp from the leading edge. Thus, from equations ( 4 ) and 

( 12 ), 


q"(x) = GcpSt^AtQ + 1 . 134Gc p Stf]Tmx (l3) 

which may be written as 

St _ 1.154 + At 0 /mx 

Strp 1 + Ato/mx ^ 

Note that, as AtQ/mx approaches zero, the Stanton number for the com- 
bined step-ramp temperature distribution approaches that of the simple 
ramp (eq. (l2a)), and as At 0 /mx becomes very large, the Stanton number 

approaches that of the isothermal plate. This is the expected limiting 
behavior . 

The step- ramp example illustrates several interesting features of 
nonisothermal heat transfer that were mentioned earlier. For example, 
if the temperature difference increases with x(m > 0), the heat-transfer 
rates will be higher than the isothermal rates. On the other hand, a 
decreasing temperature difference (m < 0) leads to lower heat transfer; 
in fact, zero heat transfer may be obtained at some point with a finite 
temperature difference, and finite heat transfer at some point with a 
zero temperature difference. Physically, these situations arise because 
the energy and consequently the temperature in the boundary layer depend 
strongly upon the ’’history” of the layer. For example, if the 
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temperature difference decreases (the plate being heated), the fluid 
near the vail is hotter than it vould have been if the plate temperature 
had been constant at its local value from the start. This means that 
the temperature gradient at the vail is less, and consequently the heat 
transfer is lover, than the isothermal values . 

Delayed ramp vail temperature . - Another case of special interest 
is that vhere the temperature difference is zero over part of the plate 
and then varies linearly, as shovn in sketch (i): 



The temperature difference is given by 


At = 0 


At = m(x - a) 
Again using equation (3), 


x < a 
x > a 


q"(x) = GcpSt 



1 



(15) 


Letting z = 1 - (£/x) 9 / 10 , the integral of equation (l5) reduces to an 
incomplete beta function. This leads to the folloving result: 


q"(x) = ^ Gc p St T (x)mxB r (8/9,10/9) (l6) 

where r = 1 - (a/x)^Al This result is of utility in approximate so- 
lutions of variable -temperature problems, as is shovn later. 


Constant temperature folloved by adiabatic vail . - Consider the 
case vhere the front part of the plate is held at constant temperature 
and the remainder of the plate is insulated, so that the heat transfer 
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is zero. The wall temperature is shown in sketch (j): 



This case is of interest in the de-icing of aircraft wings , where the 
leading edge can he heated and the warm boundary layer used to melt the 
ice from the aft portion of the wing. Similar techniques would be use- 
ful in the cooling of high-speed aircraft and missiles. 


This problem must be treated as a "specified heat input" problem. 
The heat flux over the isothermal portion can be computed from the defi- 
nition Of StrpI 


l f, (C) = Gc p St T (£)At Q x < a (17) 

The heat input is zero for x > a. Then, using the fact that St T (£) 

varies as £“ 0 - 2 , the temperature of the adiabatic wall may be determined 
from equation (6) as 


At (x) 



x > a 


(18) 

By letting z = (£/x) 9 / 10 , the integral of (l8) can be reduced to an 
incomplete beta function, leading to 


At B s (8/9,l/9) 
At 0 r(8/9)r(i/9) 


(19) 


where s = (a/x) 9 / . It should be noted that there is no way to deter- 

mine the "adiabatic decay temperature" (At) from the isothermal heat- 
transfer equation, and that nonisothermal theory is essential in the 
solution of this problem. 


P ? 7 
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sketc^fk)^^ ±Uipnt ‘ ~ Consider a step heat input as indicated in 



x 


ion 


The wall- temperature distribution rate may be determined from equat 
(6; as follows: 

(9/ 10 K f 1 r /^9/l0l' 8 / 9 

At(x) = r r i/9)r(8/9)ac p st T (x) J L 1 ' © J <|) («*) 

Z/x 

By letting z = 1 - (£/x)^/ d ^, equation (20a) integrates to give 

qoB r (1/9,10/9) 


At(x) = 


where r = 1 - (z/x) 9 / 1 ^. 


r(l/9)r(8/9)Gc p St T (x) 


(20b) 


„ 2 . ^ da y ed ^ am P heat input . - Consider a flat plate subjected to a 

delayed ramp heat input, as shown in sketch (z): 



The temperature distribution for this specif ied- heat -input problem can 
be determined using equation (6), which leads to 
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■ r ( l/9 irfe p s tl / t - (f) S/ 1' 8/9 | - D*(l) 

cl/X 


which, by appropriate substitution, reduces to 

At ^ = r(l/9)r(.8/9)Gc p St T [ B r(l/9,2°/9) - ^B r (l/9,10/9) 


( 21 ) 


where r = 1 - (a/x) 9 / 10 . This result is also of interest in the approxi- 
mate treatment of variable-heat-flux problems, as is described later. 


Approximate Methods for Handling Problems of 
Variable Wall Temperature and Heat Flux 

In many cases the prescribed wall-temperature distribution or heat 
flux will be such that the integrals resulting from the superpos it ion 
methods cannot readily be evaluated in closed form. In such cases a 
long and tedious numerical integration is required for each point at 
which the heat transfer is to be calculated. It is evident that a suit- 
able approximate technique for rapid calculations would be extremely 
useful. One possible method is to superimpose a finite number of steps 
to approximate the wall-temperature or heat-flux distribution. However, 
this leads to infinite heat-transfer rates at the discontinuities. A 
more satisfactory method is to superimpose a number of "ramps" to approx- 
imate the temperature or heat-flux distribution. Generally, a relatively 
few ramps can be used to obtain an excellent approximation to any type 
of prescribed distribution. Furthermore, since no discontinuities in 
temperature occur, no infinities in the heat -transfer rates will be ob- 
tained when ramps are used. The method of superposition of a number of 
delayed ramps is indicated by figure 3. Note that each ramp extends in- 
definitely downstream from the point of its origin. Any discontinuities 
in the prescribed distribution can be accounted for by adding a step, as 
is indicated by figure 3. Numerical examples of this method are worked 
out and compared with experiment later. 

The method indicated is similar to that described recently in ref- 
erence 13, in which some integrals of interest in laminar and turbulent 
variable- temperature problems are computed and a concise computation pro- 
cedure is presented. Computations are given for the result of both Rubes in 
(ref. 7) and Seban (ref. 14), the latter of which is similar to the result 
of the present investigation. The present scheme has the advantage that 
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the computing equations have been put in a form where loss of significant 
figures due to subtraction of numbers of the same order of magnitude is 
minimized. Moreover, the present method allows steps to be used in addi- 
tion to ramps, while in the method of reference 13 the distribution must 
be approximated entirely by ramps. Moreover, reference 13 treats only 
problems of prescribed wall temperature, while the present work allows 
calculation of both prescribed temperature-distribution and prescribed 
heat -flux problems. 

The approximation of the temperature or heat -flux distribution by 
ramps and steps may be represented mathematically by the expression 

At(x) ) * J 

| = 22 “ak ' a n) + Z# b j (22) 

q n=l j=l 

Here m n is the slope of the n th ramp and has the dimensions of °F/ft 
(or Btu/(hr)(cu ft)), and a n is its starting point; bj is the height 
of the step and has the dimensions °F (or Btu/(hr)(sq ft)); N and 

J are the number of ramps and steps, respectively, which begin upstream 
of the given location x. Differentiation of equation (22) shows that 
the slope of the approximate distribution at any point is the sum of the 
slopes of all ramps starting upstream of that point. Thus, if M(x) is 
the slope of the curve of temperature against distance at the point x, 

N 

H (x) = M(x) =^ 1 ^ (23) 

n=l 


Equation (23) allows calculation of the slopes of the component ramps. 
The step height bj is simply the rise across the step and may be cal- 
culated by a single subtraction. Thus, once the ’’break points” have 
been selected, ramps and steps may be drawn in and their parameters 
evaluated by a small amount of arithmetic. 


Prescribed temperature-difference problems . - If the temperature 
difference is represented as the sum of a number of ramps and steps, 
the heat -transfer rate is simply the sum of the heat-transfer rates due 
to the component ramps and steps. Therefore, from equations (l6) and 

(b. 


N 

q"(x) = (lO/9)Gc p St T (x)x ^ %B r (8/9,10/9) 

n=l 
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where Zj is the location of the j th step, and r n = 1 - (a n /x) 9 / 10 . 

For calculation purposes, it is convenient to rewrite equation ( 24 ) in a 
different form. In order to do this, two new functions are defined: 


and 


a (| j = ( 10/9 )B r ( 8/9, 10/9) . (1 . h 



(25) 


(26) 


Then, by using equation ( 22 ), equation ( 24 ) becomes 


q"(x) = GCpStrp(x) 


At(x) + 

_ 


X 



( 27 ) 


This . form of the computing equation has the advantage that the loss of 
significant figures due to subtraction of numbers of the same order of 
magnitude is relatively small. It also allows direct comparison with 
what would be predicted from the isothermal-plate relation (eq. ( 2 )), 
namely 3 


Qjjt = Gc p St T (x)At(x) (28) 

The effect of the nonisothermality is then concentrated entirely in the 
summation terms. In order to simplify the form, an ’’effective tempera- 
ture difference’’ At* can be defined as 


At*(x) = At(x) -h x ^ m n A ^ 


(29) 


This allows equation (27) to be written as 

q"(x) = Gc p St T (x)At* 


(30) 


The relations developed allow rapid calculation of heat-transfer 
rates in prescribed wall-temperature-difference problems. The computing 
functions, A and S, are tabulated in table III and shown in figure 4. 
An example of calculation of the heat-transfer rate for a prescribed 
temperature distribution is presented later. 


Prescribed heat-flux problems . - Prescribed heat -flux problems may 
be handled in the same manner as prescribed temperature-difference 
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problems , making use of the solution for a delayed ramp heat flux (22) 
Denoting q£(x) as the heat flux due to the n th ramp and q'!(x) as 

that due to the step heat input, equation (6) may be written as 


/ £=x 

g(£,x) 

C,=0 


N J 

E + E 


_n=0 j-0 

Exchanging the order of integration and summation gives 

•4=x ^£= x 




(31) 


Jl, /-5=x J 

At(x) = X) / A 8(«,*)<(€)« + E) 7 g(^x)q"(c)d€ (32) 

n=0 ^ =0 j=o £=° J 


But the integrals of (32) are merely those for single ramps and steps* 
therefore, from (20) and (22), 


At(x) = 


1 


r(l/9)r(8/9)Ge D St 





Br n (1/9, 20/9) 


fa. 


- lTi B r n (l/9A0/9) 


+ 2 BjV (1/9, 10/9) 


j=0 


(33) 

where r n = 1 - (a^/x) 9 / 10 , and rj = 1 - (lj/x) 9 / 10. Equation (33) may 

be put in a more useful form for computation by introducing two new 
f unct ions : 


Dl 


l-f)- 


B r (1/9, 20/9) - (f)B r (l/9,10/9) 
r(i/9)rcs/9) 


(34) 


and 


v ( i \ _ B r ( 1/9, 10/9) 

W " " ru/9)r(e/9) 


(35) 


Ihen^ using equation ( 22 ) ; equation (33) may be written 


as 


At(x) eTst^J 


(36) 


where 
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q* = q"(x) - x £ V>(|) - £ b j H (r) ( 37 ) 

n “l j=l 

and q* represents an "effective heat-transfer rate/' in which the ef- 
fects of the nonisothermality are all concentrated in the summation 
terms. Equations (36) and (37) may be used to calculate temperature dif- 
ferences in prescribed heat-flux problems. A numerical example of the 
use of equation (37) is included later with experimental verification. 

The functions D and H are tabulated in table IV and shown by figure 5. 


RESULTS AND DISCUSSION 
Comparison of Theory and Experiment 

Experiments have been performed for several cases of variable wall 
temperature and heat flux. The test apparatus consisted of a large 
heated plate with an active flow length of about 5 feet built up with 24 
individually heated strips. Within control limitations, any desired wall- 
temperature or heat-flux distribution could be obtained. The plate was 
tested in the 7-foot-diameter free- jet Guggenheim wind tunnel at Stanford 
University. Air velocities up to 130 feet per second were employed. 

This apparatus is described in detail in reference 3. 

The data from these tests have been compared with predictions ob- 
tained from nonisothermal theory. In all of these predictions, the iso- 
thermal Stanton number Stq> was computed from, the relation 

/T \~0.4 

St T Pr°- 4 = 0.0296Rex (38) 

In using this relation, all fluid properties are evaluated at the free- 
stream static temperature, and the factor (T^T^)” 0 * 4 is a temperature- 
dependent fluid-properties correction. For comparison, predictions of 
the behavior of the nonisothermal surfaces have been made employing both 
the step-function analysis presented in reference 4 and an earlier anal- 
ysis due to Rubesin (ref. 7). 

Constant heat input . - Four test runs were made in which the heat 
flux was held constant. These data are tabulated in table V(a) and are 
shown on figure 6. The data are compared with the isothermal-plate cor- 
relation (eq. (38)). the constant-heat- input predictions of the present 
analysis (eq. (lOb)), and the constant -heat -input prediction of Rubesin 
(ref. 7). Since the Stanton number for a surface with constant heat 
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input is only about 4 percent greater than that for constant wall tem- 
perature, the experimental uncertainty largely masks the predicted dif- 
ference. At the start of the plate the data appear to be closer to the 
constant -temperature correlation; this is because the strips are long 
relative to the flow length, and each strip is essentially at constant 
temperature. Farther downstream, where the strip length is small com- 
pared with the flow length, the data agree very well with the constant- 
heat- input prediction, 


St H Pr 0 - 4 = 0.0307Rex°' 2 &) (39) 

' 00 / 

The prediction of Rubesin (Stjj/St^ = 1.06) appears to be slightly high 
over the entire test range. 

Doubl e- step wall temperature . - One test run was made for which the 
wall- temperature distribution was a "double step, " a discontinuity in the 
wall temperature occurring in the center of the plate. These data are 
tabulated in table V(b) and shown in figure 7. Over the first portion of 
the plate, where the plate is at constant temperature, the predictions of 
the present analysis and of Rubesin (ref. 7) are identical. After the 
discontinuity, Rubesin's analysis is high, while the prediction of the 
present analysis is quite good. The predicted heat- transfer rates were 
determined by adding the rates due to the step at the leading edge to 
those due to the step downstream. 

Step-ramp wall temperature . - One test run was made in which the 
wall temperature was varied linearly from the leading edge. Because of 
control limitations, it was necessary to have a small step in the wall 
temperature at the leading edge . The data from this run are tabulated 
in table V(c) and shown in figure 8. The linearity of the wall tempera- 
ture was quite good. The heat-transfer data are compared with the pre- 
dictions from the present analysis, with the result of Rubesin (ref. 7), 
and with what would be obtained from use of the uncorrected isothermal- * 
plate correlation (eq. ( 38)), employing the local temperature difference. 
The predictions of the present analysis, which were determined from, equa- 
tion (13), agree very well with the experimental values. The heat-flux 
predictions based on the analysis of Rubesin are high, while the iso- 
thermal relation (38) predicts heat fluxes that are slightly low. How- 
ever, the ramp obtainable with the experimental apparatus was not very 
steep, and thus the effect of the nonisothermality is not as great as 
might be obtained (see eq . (14)). In fact, use of the isothermal equa- 
tion (38) and the local temperature differences gives results that may 
well be satisfactory as first approximations for engineering design 
calculations . 

Constant temperature followed by adiabatic wall . - One test run was 
made in which the forward part of the plate was held at constant tempera- 
ture and the power was turned off on the last 12 strips. The data from 
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this run are tabulated in table V(d) and shown by figure 9. In this ex- 
ample, the hot air flowing over the unheated portion of the plate warms 
the plate, and the result of interest is the wall-temperature distribu- 
tion downstream of the discontinuity in heat flux. This situation is of 
interest in the de-icing of wings, where the leading edge may be heated 
and the hot air used to melt the ice from the after portion of the wing. 
Similar techniques are of interest in the cooling of missiles, where the 
front of the missile is cooled and the cold boundary layer is used to 
cool the rearward sections. The data are compared with the temperatures 
predicted by the present analysis (eq. (l9)) and a similar equation ob- 
tained from the Rubesin analysis (ref. 7). Again the Rubesin prediction 
is high, while the agreement between the data and the present analysis 
is quite good. The measured temperatures are expected to lie slightly 
higher than the predictions because of conduction in the test plate. 

This conduction effect is greatest just downstream, of the discontinuity 
in heat flux, and agreement of the data with the predicted temperature 
distribution is poorest in this region. The over-all agreement of the 
data with predictions is considered quite satisfactory. 

Double -pulse heat input . - One test run was made with a double or 
step heat-flux distribution. The data from this run are tabulated in 
table V(e ) and shown by figure 10. Because of conduction in the plate, 
it was not possible to obtain sharp discontinuities in the heat flux, as 
the data Indicate. In making the predictions for the wall temperature, 
the actual nonperfect pulses were approximated by perfect pulses, so that 
the total heat transfer was the same. The approximate pulses are indi- 
cated in figure 10. 

The predicted temperatures were determined by superposition of three 
steps of heat input, using the step heat-input results (eq. (20)). The 
agreement between the predicted and experimental temperatures is excel- 
lent, except near the ends of the pulses. The over-all agreement is con- 
sidered quite satisfactory. 

Irregular wall temperature . - One test run was made in which an ex- 
tremely irregular wall-temperature distribution was maintained in the 
plate. The data from this run are tabulated in table V(f) and shown in 
figure u( a) . 

The experimental heat-transfer rates are compared with heat-transfer 
rates predicted from the approximate methods described earlier. Figure 
11 (b) shows how the wall-temperature distribution was approximated by 
seven ramps and one step. Again the predicted heat-transfer rates are 
in excellent agreement with the predictions of the approximate method. 

The data for strip 2 are low, probably because the flow was not fully 
turbulent over this strip. The heat-transfer rate predicted by use of 
the isothermal-plate relation (38) is in considerable error, especially 
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where the wall temperature is changing rapidly. However, the values are 
of the same order of magnitude, and it can therefore be concluded that 
the isothermal equation may be used to obtain a first estimate of the 
heat transfer from nonisothermal surfaces . 

This same example has also been treated as a "prescribed-heat-flux" 
problem, using the approximate methods described earlier. Figure 12(a) 
shows how the heat flux was approximated by five steps and six ramps. 
Figure 12(b) shows the temperature differences computed by the approxi- 
mate methods compared with the experimental temperatures . The agreement 
near the leading edge is not too good, but this can be attributed to the 
approximations made on the heat -flux distribution, and to the fact that 
here the flow may not have been fully turbulent . Farther downstream, 
the predicted and experimental temperatures agree very well. The pre- 
dictions made from the isothermal equation (38) are in considerable er- 
ror, especially where the heat flux is changing rapidly. They are, how- 
ever, of the same order of magnitude, and thus it appears that the iso- 
thermal equation can also be used to obtain a first approximation to the 
temperature difference in a prescribed heat-flux problem. 


Determination of Heat-Transfer Rates by Approximate Methods 

In order to illustrate the approximate methods for handling pre- 
scribed wall-temperature problems, the "irregular" example described in 
the preceding section will be worked out in detail. The test data for 
this run are shown in figure 11. The prescribed temperature difference 
will be approximated by seven ramps and one Step, as is indicated by 
figure ll(a) . 

The first step is to determine the location and slopes of the ramps. 
The "break points" for the ramps and the temperature difference for the 
step will be taken as indicated in the following table, where n denotes 
the number of the ramp starting at a n : 


n 

a n , 

ft 

At(a n ) , 
°F 

1 

0 

4.0 

2 

.70 

20.2 

3 

1.20 

11.4 

4 

1.82 

7.6 

5 

3.10 

18.4 

6 

3.50 

24.6 

7 

4.00 

16.7 

end 

4.80 

13.2 
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The slope of each component ramp may he determined readily , since 
slope of the approximate temperature distribution at any point is 
the sum of the slopes of all ramps starting upstream. Therefore 
the first ramp, 


ra l 


20.2 - 4.0 
0.70 - 0 


= 23.14 °F/ft 


To find the slope of the second ramp, 

m i + m 2 = Mrffi s = - 17 - 60 ° F / ft 

m 2 - -17.60 - 23.14 = -40.74 °F/ft 
Similarly, for the third ramp, 


m l + m 2 + m 3 


7.6 - 11.4 
1.82 - 1.20 


-6.13 °F/ft 


m 3 = -6.13 + 40.74 - 23.14 = 11.47 °F/ft 


For the fourth ramp, 


rri4 + (-6.13) 


18.4 - 7.6 
3.10 - 1.82 


8.44 0 F/ ft 


m 4 = 14.57 °F/ft 


For the fifth ramp, 


^ + 8 - 44 = M-tfm = 15.50 o F/ft 

m 5 - 7.06 °F/ft 

For the sixth ramp, 

“6 + 15 - 50 = i!oO - 5 4 50 “ - 15 - 80 ° F / ft 
m g = 31.30 °F/ft 
Finally, for the seventh ramp, 


the 

simply 

for 
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m y + (-15.80) = g g§-5 ^ = -4.37 <*/ft 
m 7 = 11.43 °F/ft 

* « 

In addition to the seven ramps , the superposition involves a step of 
4.0° F at the leading edge, where = 0. Summarizing, 


| Ramps 

n 

a n.> 

m n> 


ft 

°F/ft 

1 

0 

23.14 

2 

.70 

-40.74 

3 

1.20 

11.47 

4 

1.82 

14.57 

5 

3.10 

7 .06 

6 

3.50 

-31.30 

7 

4.00 

11.43 


Steps 

j 

ft 

b j’ 
°F 

1 

0 

4.0 


At this point a check should he made. For 4.00 < x < 4.80 feet 
the temperature distribution should be given by 

At(x) = 4.0 + 23.14(x - 0) - 40.74(x - 0.70) + 11.47(x - 1.20) 

+ 14.57 (x - 1.82) + 7.06(x - 3.10) - 31.30(x - 3.50) + 11.43(x - 4.00) 

By substituting x - 4.80 feet, the temperature at the "end" may be cal- 
culated from this relation as 

At (4 . 80 ) = 4.0 + 111.07 - 166.93 + 41.29 + 43.42 + 12.00 

-40.69 + 9.14 = 13.3° F 


This agrees with the value that exists at x = 4.80 feet, and therefore 
the calculated ramp slopes close with the proper value, providing a 
check . 


The calculation of the heat -transfer rates can best be handled by 
tabular computation. For example, the calculation of the heat-transfer 
rate at x = 3.5 feet is as follows: 

First, tabulate the parameters of interest. From the data, 
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x = 3.5 ft 

G - 32,400 lb/(hr)(sq ft) 
[ i = 0.0439 lb/ (hr) (ft) 
c p = 0.24 Btu/(lb)(°F) 


At - 24.6° F 
Pr - 0.70 
T^ = 528° R 
T w = 553° R 


Next, in a tabular manner, calculate the sums appearing in equation (27): 


n 

% 

a n /x 

A(a n /x) 

“n 

A(a n /x)m n 

1 

0 

0 

0.135 

23.14 

3.12 

2 

.70 

.200 

.130 

-40.74 

-5.30 

3 

1.20 

.343 

.124 

11.47 

1.42 

4 

1.82 

.520 

.112 

14.57 

1.63 



.886 

.053 

7.06 

.37 





^A(a n /x.)m n = 1.24 °F/ft 

j 



S(Zj/x) 


S(7j/x)bj 

1 

0 

0 

0 

4.0 

0 





£s(lj/x)bj = o° F 


Note that only the ramps and steps starting upstream of x are used in 
these calculations. The functions A and S are determined from fig- 
ure 4, which may be replotted from the values in table III. The '’effec- 
tive" temperature in equation (29) is therefore 

At* = At(x) + x^m n A(a n /x) + 10 j s (7j/x) 

- 24.6 + 3.5(1.24) + 0 

= 28.9° F 

The Reynolds number is 

Re x = Gx/jjl - (32,300)(3.5)/0.0439 
= 2.58X10 6 

The isothermal Stanton number Stiji may be determined from equation (38) 
as 
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St T = 0. 0296(2. 58X10 6 )" 0,2 (0.7)"‘ 0 * 4: (553/528)"‘ 0 - 4 
= 0.00176 

Finally, the heat -transfer rate may he calculated from equation (30 ) as 
q" (3.5) = (32,300) (0.24) (0 .00176 ) (28 . 9 ) 

= 394 Btu/(hr)(sq ft) 

Similar calculations may he repeated at any point for which the heat- 
transfer rate is to he predicted. In making the predictions for this 
example, heat-transfer rates were calculated at the six "break points," 
the "end" (4.80 ft), and at x = 0.4 foot. The total calculating time 
was slightly less than 1 hour. 


Determination of Wall Temperatures hy Approximate Methods 

In order to illustrate the approximate methods described earlier 
for handling prescribed heat-flux problems, the "irregular" example de- 
scribed earlier will be worked out in detail, assuming the heat flux is 
known and the wall temperature is to be calculated. The data from this 
run are shown in figure 12 . The prescribed heat flux will be approxi- 
mated by six ramps and five steps, as is indicated by figure 12(a). It 
will be assumed that the heat -transfer rate was constant from, the lead- 
ing edge, having a value equal to that measured for the second strip. 
This assumption will introduce some error at the start, but the boundary 
layer on the first strip is probably laminar or transitional, and thus a 
more accurate calculation is not possible. 


The first step is to evaluate the parameters in the step-ramp ap- 
proximation. The heat-flux discontinuities will occur midway between 
data points; at these points the two adjacent strips are separated by a 
thin insulator. The discontinuity and break-point locations are as 
follows : 


ft 

q"(x), 

Btu/(hr)[sq ft) 

0 - 

0 

0+ 

162.5 

.46- 

162.5 

.46+ 

372.5 

.90- 

372.5 

.90+ 

235.0 

1.10- 

235.0 

1.10+ 

175.0 

1.70 

110.0 

2.16- 

110.0 

2.16+ 

152.5 

3.00 

280.0 

3.50 

400.0 

3.88 

245.0 

4.80 

130.0 
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First, the height of each step is calculated as follows: 

hi = 162.5 - 0 = 162.5 Btu/(hr)(sq ft) Zq = 0 ft 

h 2 = 372.5- 162.5 = 210.0 Btu/(hr)(sq ft) Z 2 = 0.46 ft 

h 3 = 235.0 - 372.5 = -137.5 Btu/(hr)(sq ft) Z^ = 0.90 ft 

hq = 175.0 - 235.0 = -60.0 Btu/(hr)(sq ft) Z^ = 1.10 ft 

b 5 = 152.5 - 110.0 = 42.5 Btu/(hr)(sq ft) Zg = 2.16 ft 

The slope of each ramp is calculated from equation (23): 

m x = (110.0 - 175 .0)/(l . 70 - 1.10) 

= -108.3 Btu/ (hr) (cu ft) 

m 2 + (-108.3) - ( 110.0 - 110.0)/(2.16 - 1.70) = 0 
m 2 = +108.3 Btu/(hr)(cu ft) 
m 3 + 0 = (280.0 - 152.5)/(3.00 - 2.16) - 151.8 
m 3 = 151.8 Btu/(hr)(cu ft) 
m 4 + 151.8 = (400.0 - 280.0)/(3.50 - 3.00) = 240.0 
m 4 = 8Q-2 Btu/(hr)(cu ft) 

m 5 + 240.0 = (245.0 - 400.0)/(3.88 - 3.50) = -407.9 
m 3 = -647.9 Btu/(hr)(cu ft) 

m 6 + (-407.9) = (130.0 - 245.0)/(4.80 - 3.88) = -125.0 
m^ = 282.9 Btu/(hr)(cu ft) 
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The foregoing calculations may be summarized as follows: 


Ramps 

n 

a n> 

m n> 


ft 

Btu/(hr)(cu ft) 

1 

1.10 

-108.3 

2 

1.70 

108.3 

3 

2.16 

151.8 

4 

3.00 

88.2 

5 

3.50 

-647.9 

6 

3.88 

282.9 


Steps 

j 

V 

b j> ^ 


ft 

Btu/(hr)(sq ft) 

1 

0 

162.5 

2 

.46 

210.0 

3 

.90 

-137.5 

4 

1.10 

-60.0 

5 

2.16 

42.5 


At this point, it is desirable to make a check. At x = 4.80 feet, the 
heat flux is given by 

q"(x) = 162.5 + 210.0 - 137.5 - 60.0 + 42.5 - 108.3(4.80 - 1.10) 

+ 108.3(4.80 - 1.70) + 151.8(4.80 - 2.16) 

+ 88.2(4.80 - 3.00) - 647.9(4.80 - 3.50) + 282.9(4.80 - 3.88) 

= 130.0 Btu/(hr)(sq ft) 

This result agrees with the value of q" at the ’’end, " and thus it ap- 
pears that no errors have been made in determining the b and m. 

The calculation for the temperature difference at various points on 
the plate can be handled in a tabular manner. To illustrate the method, 
the calculation for the temperature difference at x = 3.5 feet follows. 

First, tabulate the parameters of interest. From the experimental 

data, 


x = 3.5 ft 

G = 32,400 lb/(hr)(sq ft) 
[i = 0.0439 lb/(hr)(ft) 


c p = 0.24 Btu/(lb)(°F) 
Pr - 0.70 

q" = 400 Btu/(hr)(sq ft) 


Next, in a tabular manner, calculate the sums appearing in equation (37): 
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Steps 


j 

h 

Jj/x 

H(Zj/x) 

h 

Y ( V X > 

1 

0 

0 

0.041 

162.5 

6.66 

2 

.46 

.131 

.055 

210.0 

11.55 

3 

.90 

.257 

.070 

-137.5 

-9.63 

4 

1.10 

.314 

.080 

-60.0 

-4.80 

5 

2.16 

.617 

.133 

42.5 

5.65 





^bjHpj/x) = 9.43, 





Btu/(hr)(sq ft) 


Ramps 

n 

a n 

a n /x 

D ( a n / X ) 

m n 

% D (an/ x ) 

1 

1.10 

0.314 

0.116 

-108.3 

-12.56 

2 

1.70 

.486 

.100 

108.3 

10.83 

3 

2.16 

.617 

.085 

151.8 

12.90 

4 

3.00 

.857 

.043 

88.2 

3.79 





EvXvx) = 14 - 96 




x fX D K/ x ) 

= 52.36 Btu/(hr)(sq ft) 


Therefore, the effective heat flax q* is 

q* - 400.0 - 52.36 - 9.43 
= 338.2 Btu/(hr)(sq ft) 


The Reynolds number is 


Re x = Gx/|i = (32,300)(3.5)/0.0439 
= 2.58X10 6 

The isothermal Stanton number St T may be determined from equation (38). 
However, since the temperature is not known, the fluid-properties correc- 
tion factor (T^/t^) cannot be evaluated. This is, however, a small 

correction, which may be neglected without serious error. Thus, 

St T - 0. 0296(2. 58Xl0 6 )-°* 2 (0.7)-°- 4 


= 0.00179 
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Finally, the temperature difference may be calculated from equation (36): 
At (3.5) = 338.2/(32, 300)(0.24)(0. 00179) 

= 24.4° F 

Similar calculations may be repeated at each point for which the tempera- 
ture difference is desired. One must be careful, however, at the steps. 
Just upstream of a step (x = Z“j), the step should not be included in the 

calculation. But just downstream of the step (x = zip, the step must be 

included, or else the computations would indicate a discontinuity in the 
wall temperature. Bote that use of the isothermal equation (38) to pre- 
dict the temperature leads to considerable error in heat flux at the 
steps where this relation indicated a discontinuity in the wall tempera- 
ture. However, the isothermal equation is useful for an order-of- 
magnitude check on. the calculations. 


SUMMARY OF RESULTS 

? The good agreement between predicted and measured values of both 

heat -transfer rates and temperature distributions for the variety of 
situations shown is believed to substantiate the present theory for non- 
isothermal surfaces entirely adequately. Consequently, theory should be 
sufficiently accurate to predict heat transfer from nonisothermal sur- 
faces in virtually any similar type of situation. 

The results also show that the usual correlation for the heat- 
transfer rates, that in the solution for the isothermal surface, is in 
fair agreement with the data in all cases except where rapid changes in 
wall temperature or heat flux occur. Hence, in many cases the simpler 
isothermal equations can be used to obtain a first approximation for the 
heat-transfer rates (or temperature differences), even though the surface 
is not isothermal. 

In cases where rapid changes in wall temperature occur, where a 
fraction of the surface is adiabatic, or where high accuracy is required, 
the nonisothermal theory should be used. In these situations, the ap- 
proximate methods presented herein provide a simple, rapid computation 
method of good accuracy for the prediction of the temperature distribu- 
tion or heat-transfer rates for any arbitrarily prescribed conditions. 


Stanford University, 

Stanford, Calif., October 22, 1957. 
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TABLE I. - INCOMPLETE BETA FUNCTIONS FOR TURBULENT 


NONISOTHERMAL HEAT TRANSFER 


r 

B r (8/9, 1/9) 

B r (8/9, 10/9) 

B r (l/9,10/9) 

B r (l/9,20/9) 

0 

0 

0 

0 

0 

.1 

.1518 

.1445 

6.9604 

6.8837 

.2 

.2948 

.2661 

7.5087 

7.3445 

.3 

.4457 

.3792 

7.8447 

7.5898 

.4 

.6109 

.4863 

8 . 0886 

7.7414 

.5 

.7971 

.5886 

8.2796 

7.8401 

.6 

1.0151 

.6864 

8.4357 

7 • 9046 

.7 

1.2833 

.7797 

8.5664 

7.9456 

• 8 

1.6412 

.8682 

8.6772 

7.9694 

.9 

2.2089 

.9505 

8.7702 

7.9809 

1.0 

9.1853 

1.0206 

8.8439 

7.9839 

r 

B r (l/9,8/9) 

B r ( 1/9,10/9) 

B r (l/9,20/9) 


0 

0 

0 

0 


.01 

5.3960 

5.3948 

5.3887 


.02 

5.8286 

5.8260 

5.8131 


.03 

6.0979 

6.0938 

6.0735 


.04 

6.2967 

6.2910 

6.2631 


.05 

6.4555 

6.4482 

6.4125 


.06 

6.5883 

6.5794 

6.5357 


.07 

6.7029 

6.6923 

6.6406 


.08 

6.8039 

6.7916 

6.7316 


.09 

6.8943 

6.8802 

6.8119 


.10 

6.9764 

6.9604 

6.8837 












mBLE II. - SUMMARY OF TURBULENT NONISOTHERMAL HEAT- TRANSFER 

SOLUTIONS FOR A FLAT PLATE 


SPECIFICATION 


ARBITRARY TEMPERATURE 


ARBITRARY HEAT FLUX 


TEMPERATURE 


At = A+( x ) 


p T y-0 


CONSTANT TEMPERATURE 
At 

n — 


STEP TEMPERATURE 
At 1 

irTb 

o 1 1 — 1 


RAMP TEMPERATURE 


At = 0 X£ l 

At = b x>Z 


At = 0 * <0 

At = m(x -a) x>a 


CONSTANT 

ADIABATIC 

TEMPERATURE- 

WALL 

At = b 

X * 

At h i\ 
n b a ' 

77777T777i 

r+" X 

At * 0.1084 b B s (8/9,l/9) 

X2 


CONSTANT HEAT INPUT 


STEP HEAT INPUT 


RAMP HEAT INPUT 



’ GcpSt T 


At =0 

0.959 b B r (l/9,IO/9) 
GcpSt T Bi (l/9 ,10/9 ) 


At = 0 

«+ 0. 1084 mx f D.fi/Q _ In 


X<1 

X£Z 


X<a 


fltiLi 084 .™.* f B r (1/9,20/9) -(d/x)B,(l/9, 10/9)1 

x Gc p St T L X S Q J 


HEAT FLUX 


q" = Gc P St T J(l-(^/x) 9/l °] l/9 dAt(p 
bo 


q"=q"(x) 


q" 5 Gc P St T b 


q" = 0 X<1 

v- 1/9 

q”= Gc p St T b[i- (Z/x) 9/l0 I Xfcl 


q" = 0 x<a 

q" = unmxGcpSt-j- B r (e/9,i0/9) x>a 


q ,,= G CpSt T b X<a 


q" -0 X < l 

q” = b x a. I 


II ~ 

q = 0 x<0 

q" = m(x-a) X>a 


St T = 0.0296 Pr -0 ’ 4 Re x 02 


r * !'>!'/>;« 

- 1 - U/x) 


s = (a/x) 



























a/x 

A (a/x) 

l/x 

s (l/x) 

0 

0.135 

0 

0 

.1 

.133 

.1 

.015 

.2 

.130 

.2 

.030 

.3 

.126 

.3 

.047 

.4 

.121 

.4 

• 066 

.5 

.113 

.5 

.089 

.6 

.103 

.6 

.116 

.7 

.100 

.7 

ol52 

.8 

o074 

.8 

.208 

.9 

.050 

.9 

.300 

1.0 

0 

1.0 

oo 


TABLE IV, - FUNCTIONS FOR PRESCRIBED 


HEAT-FLUX CALCULATIONS 


a/x 


.1 

.2 

.3 

.4 

0 5 

.6 

.7 

.8 

,9 

1.0 


D(a/x) 



H (l/x) 


0.135 

.130 

.124 

.117 

.109 

.099 

.087 

.072 

.056 

.033 


.1 
.2 
.3 
.4 
.5 
. 6 
.7 
.8 
.9 
1.0 


0.041 
.050 
.056 
.078 
.092 
.110 
.130 
.158 
.194 
„ 256 
1.000 


0 





































TABLE V. - Concluded. EXPERIMENTAL HEAT-TRANSPER DATA 










WEEK 


LOAD DISTRIBUTION 
w(x) V 


T w (x) 


WALL - TEMPERATURE 
DISTRIBUTION, T w (x)— * 



<p - f W (x)<d ( 1, x) dx 

CJ (1, x) = DEFLECTION AT 7 
DUE TO UNIT LOAD AT x 


q"(x) 


q"CO = f VI, *) d V x > 

•'o 

h T (7,x) = HEAT - TRANSFER 
RATE AT 7 DUE TO UNIT 
TEMPERATURE RISE AT x 


Figure 1. - Analogy of superposition techniques for beam deflection and heat transfer 








Cif'-b' back 


m l + m 2 +m 3 


m, + m 2 


nrij +m2 + 
m ^ +■ bn ^ 


prescribed temperature 

OR HEAT FLUX 


APPROXIMATION BY STEP 
STEPS AND RAMPS ~T l 
, b. I 


° 3 ^ 




m 


Figure 3. - Approximation of prescribed temperatures or heat flux by steps and ramps 




.12 


.10 
.08 

A (a/ x) 

.06 

.04 

.02 

| 

0 .2 .4 .6 .8 

a/x 


Figure 4. - Functions 
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HEAT FLUX PREDICTIONS 

C RUBESIN — 

J PRESENT INTEGRAL 
] ANALYSIS “ 

^ ISOTHERMAL __ 

EQUATION 


< 200 

LlJ 

X 


O MEASURED TEMPERATURES 
A MEASURED HEAT FLUX 


DISTANCE FROM LEADING EDGE, x,ft 

Figure 7. - Example for double-step temperature distribution. 




HEAT FLUX, Btu / (hr)(sq ft) 
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CD 


Figure 8. - Example for step-ramp wall temperature. 








TEMPERATURE DIFFERENCE, At, 



HEAT FLUX, Btu / (hr)(sq ft) 
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Figure 11. - Prescribed temperature treatment of irregular example. 




A MEASUREO HEAT FLUX 


APPROXIMATION BY 
STEPS AND RAMPS 


O MEASURED TEMPERATURES 



T 

0 

j 

z: 


f 



PREDICTED 

TEMPERATURES 


.TEMPERATURES PREDICTED FROM 
ISOTHERMAL EQUATION 


DISTANCE FROM LEADING EDGE, X, ft 

Figure 12. - Prescribed heat-flux treatment of irregular example. 
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